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Abstract 

Relativistic jets from compact objects are ubiquitous phenomena in the Unvierse, but their driving mechanism 
has been an enigmatic issue over many decades. Two basic models have been extensively discussed: magnetohy- 
drodynamic (MHD) jets and radiation-hydrodynamic (RHD) jets. Currently, the former is more widely accepted, 
since magnetic field is expected to provide both the acceleration and collimation mechanisms, whereas radiation 
field cannot collimate outflow. Here, we propose a new type of jets, radiation-magnetohydrodynamic (RMHD) jets, 
based on our global RMHD simulation of luminous accretion flow onto a black hole shining above the Eddington 
luminosity. The RMHD jet can be accelerated up to the relativistic speed by the radiation-pressure force and is col- 
limated by the Lorentz force of a magnetic tower, inflated magnetic structure made by toroidal magnetic field lines 
accumulated around the black hole, though radiation energy greatly dominates over magnetic energy. This magnetic 
tower is collimated by a geometrically thick accretion flow supported by radiation-pressure force. This type of jet 
may explain relativistic jets from Galactic microquasars, appearing at high luminosities. 

Key words: accretion, accretion disks — black hole physics — magnetohydrodynamics: MHD — radiative 
transfer — ISM: jets and outflows 



1. Introduction 

Astrophysical jets, highly collimated, high velocity plasma 
outflows, are found in many astrophysical objects, including 
young stellar objects (YSOs), microquasars, active galactic nu- 
clei (AGNs), and gamma-ray bursts (GRBs). Importantly, the 
jets (except for YSO jets) are, at least, mildly relativistic and 
sometimes highly relativistic and are collimated with an small 
opening angle which is only a few degrees. It is widely ac- 
cepted from the observations that these relativistic jets are as- 
sociated with accretion flows onto a central compact object, 
such as black-holes (see review by Mirabel & Rodriguez 1999). 
Despite a long research history since the 1970's, how to col- 
limate and accelerate astrophysical jets still remains an open 
question (e.g., Rees 1984; Meier 2003). 

There are two models proposed and considered for the 
driving mechanism of relativistic jets: magnetohydrody- 
namic (MHD) jets which are driven by the magnetic process 
(Lovelace 1976; Blandford & Payne 1982; Uchida & Shibata 
1985; Hawley & Balbus 2002) and radiation-hydrodynamic 
(RHD) jets which are driven by the strong radiation-pressure 
force (Bisnovatyi-Kogan & Blinnikov 1977; Eggum et al. 
1988; Sikora et al. 1996; Okuda et al. 2009). The former is 
very promising, since not only the acceleration but also the col- 
limation can be accounted for, whereas the latter cannot gen- 
erally explain the collimation. However, we should remark 
that most of the previous MHD disk-jet simulations assumed 
that the underlying accretion flows are radiatively inefficient; 
i.e., low-luminosity accretion flow (Ichimaru 1977; Narayan & 
Yi 1994). By contrast, radiation field dominates the dynamics 
of luminous (radiation-dominated) accretion flows, from which 



highly relativistic jets are occasionally observed, especially in 
the cases of Galactic microquasars. 

Here, we examined the complex interactions between mat- 
ter, radiation, and magnetic fields in the high-luminosity 
accretion flow and jet, based on our global radiation- 
magnetohydrodynamic (RMHD) simulations. Our goal is to 
give good physical explanation why and how relativistic jets 
can be produced in luminous accretion flows and what kinds of 
roles magnetic field plays in the situations that radiation field 
energy dominates over matter and magnetic energy. 

2. Overview of the RMHD Simulations 

We analyzed the RMHD simulation data of Ohsuga et al. 
(2009). Here, we outline the numerical method. We extended 
the resistive MHD simulations performed by Kato et al. (2004), 
and solved the equations including the radiation transfer terms 
numerically in a 2.5-dimensional domain. We used cylindrical 
coordinates (r, 0,z), where r, 9, and z are the radial distance, 
the azimuthal angle, and the vertical distance, respectively. The 
basic equations are the continuity equation, 
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the energy equation of the radiation, 
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Here, p is the gas mass density, v = {v r ,v v ,v z ) is the flow 
velocity, c is the speed of light, e is the internal energy den- 
sity of the gas, p gas is the gas-pressure, B = (B r ,B v ,B z ) is 
the magnetic field, J{= cV x B/4ir) is the electric current, 
Bbb is the blackbody intensity, Eq is the radiation energy den- 
sity, Fq is the radiative flux, Po is the radiation-pressure ten- 
sor, k is the absorption opacity, and \ is the total opacity, re- 
spectively. We adopted anomalous resistivity, r\ & (Yokoyama & 
Shibata 1994), the flux-limited diffusion (FLD) approximation 
to evaluate Fq and Po (Levermore & Pomraning 198 1), and the 
pseudo-Newtonian potential, V'PN; to incorporate the general 
relativistic effects, given by ^pn = —GMbr/{V r 2 + z 2 — r s ) 
(Paczyhski & Wiita 1980). Here, r s (= 2G¥ BH /c 2 ) represents 
the Schwarzschild radius, where G is the gravitational constant 
and Mbh is the mass of the black hole. As the equation of state, 
we use the relation, p gas = (7— l)e = pk-QT g&s / pra p , where 
7 is the specific heat ratio, k-g, is the Boltzmann constant, T gas 
is the temperature of the gas, \i themean molecular weight, and 
m p is the proton mass, respectively. 

The computational domain is 3r s < r < 103 r s and < z < 
91.6 r s . The grid spacing of the radial distance, Ar, and the 
vertical distance, Az, is uniform (A?' = Az = 0.2 r s ). The re- 
flection symmetry relative to the equatorial plane (with z = 0) is 
assumed. Free boundary conditions for matter and a magnetic 
field are adopted; i.e., the matter can freely go out, but cannot 
to come in, and the magnetic field cannot change across the 
boundaries. Radiation is set to go out of the boundaries with 
the radiative flux of cEq, except at the innermost boundary at 
r = 3r s and z > 3r s , through which no radiation can go out. We 
also set a non-rotating stellar-mass black hole (Mbh = 10M Q ) 
at the origin, and initially put a rotating torus at radii around 
40 r s threaded with closed poloidal magnetic field (plasma- 
j3 = Pgas/Pmag = 100), and a non-rotating isothermal corona. 
Here, p mag (= B 2 /8tt) is magnetic pressure force. 

3. RMHD Jet Emerging from a Luminous Accretion Flow 

In Figure 1, we show a bird's-eye view of luminous accre- 
tion flow (at the center, with the brown color), outflow (with 
the white and the blue color), and the associated magnetic field 
lines (with the white lines). In the jet region only the parts 
with high velocities (~ 0.6c — 0.7c) are illustrated by the blue 
color. The emergence of a magnetic tower formed by the ac- 
cumulation of the toroidal magnetic field which emerged from 
the underlying accretion flow is evident. This is quite rem- 
iniscent of non-radiative MHD jets (Lynden-Bell 1996; Kato 



et al. 2004; Nakamura et al. 2006), in spite of large radiation 
energy density which by far exceeds the magnetic energy in 
the present case (discussed later). From the simulation data, 
we obtain the following quantities; the mass-outflow rate is 
Mout ~ 10Le/c 2 ; the photon luminosity is i p h ~ L E ; the ki- 
netic luminosity is Lkin ~ O.ILe; the Poynting luminosity is 
ipoy ~ O.OILe- Note that the Eddington luminosity (Le) can 
be exceeded in the case of disk accretion onto black holes be- 
cause of significant anisotropy of the radiation field and pho- 
ton trapping effects (Ohsuga & Mineshige 2007). The mass- 
accretion rate to the central black hole is M acc ~ 100Le/c 2 , 
that exceeds the critical rate giving rise to the Eddington lu- 
minosity, on the order of IOXe/c 2 , i.e., the simulated flow 
is supercritical (Shakura & Sunyaev 1973; Abramowicz et al. 
1988; see Chap. 10 of Kato et al. 2008 for a review). Here, the 
mass-accretion rate and the mass-outflow rate are calculated by 
summing up the mass passing through the inner boundary and 
the upper boundary per unit time with higher velocities than 
the escape velocity. Photon luminosity, kinetic luminosity, and 
Poynting luminosity are calculated based on the radiative flux, 
the kinetic energy, and magnetic energy at the upper boundary. 

3.1. Acceleration Mechanism of RMHD Jet 

To see what force is responsible for the acceleration of the 
jet, we evaluate the time-averaged strengths of the vertical 
components of the gravitational force, the gas-pressure force, 
the radiation-pressure force, and the Lorentz force at a radius, 
r = 10 r s , and display them in the left panel of Figure 2 as 
functions of the vertical coordinate, z. This figure unambigu- 
ously shows that the RMHD jet is accelerated by the radiation- 
pressure force (indicated by the red line), and that the Lorentz 
force (indicated by the blue line) does not contribute much to 
the acceleration. This is a natural consequence, since the radi- 
ation energy density greatly exceeds the magnetic energy den- 
sity, typically, by more than one order of magnitude. To under- 
stand the direction of the radiation-pressure force, we display 
in the left wall surface of Figure 1 the color contours of the zz- 
component of the radiation-pressure tensor, Pg 2 , overlaid with 
the radiation-pressure force vectors (arrows), xFq I P c > on the 
meridional plane. We can clearly see that the vertical, outward 
component of the radiation-pressure force dominants in the jet 
region (around the z-axis), since the radiation energy density 
rapidly varies in the vertical direction, but not so much in the 
horizontal direction; that is, the significant anisotropy of the 
radiation field (Ohsuga & Mineshige 2007). 

3.2. Collimation Mechanism of RMHD Jet 

Next, we examine the collimation mechanism of the RMHD 
jet in the same manner as we did for the acceleration mecha- 
nism. We calculate the time-averaged strengths of the radial 
components of the gravitational force, gas-pressure force, cen- 
trifugal force, radiation-pressure force, and the Lorentz force 
at z = 40 r s , and display them in the right panel of Figure 2 
as functions of r. In contrast with the left panel, the Lorentz 
force (indicated by blue solid line) dominates; i.e., the mag- 
netic force works as the a collimation force in the RMHD jet 
against the outward (centrifugal) force (indicated by the ma- 
genta dashed line). We note that magnetic -pressure force and 
the tension force are comparable. To understand the direction 



A Novel Jet Model: Magnetically Collimated, Radiation-Pressure Driven Jet 



3 




Fig. 1. Bird's-eye view of the luminous accretion flow and the associated RMHD jet. The accretion flow (the gas mass density, brown) and the RMHD 
jet in which velocities exceed the escape velocity (the velocity, white, blue) are plotted. The high-speed jet (~ 0.6c — 0.7c) is represented by blue. White 
lines indicate the magnetic field lines. The zz-component of the radiation-pressure tensor (color), Pq z . overlaid with the radiation-pressure force vectors 
(arrows) on the meridional plane is projected on the left wall surface, while the magnetic pressure from the azimuthal component of the magnetic field 
(color), B^/8n, overlaid with the Lorentz force vectors (arrows) on the meridional plane is projected on the right wall surface. The color bar corresponds 
to the logarithmic value of Pq z and B^/8ir. The arrows on each wall are displayed only in the region in which their values are larger than 10 10 dyn/g. 
Each quantity is time-averaged over 1 s, which corresponds to the accretion timescale at several tens of the Schwarzschild radius (r s ). 
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Fig. 2. Acceleration and collimation mechanism of the RMHD jet. Left panel: The vertical profiles of the gravitational force (black), the gas-pressure 
force (green), the radiation-pressure force (red), and the Lorentz force (blue) at r = 10 r E . The solid lines and the dashed lines indicate the outward force 
and the inward force, respectively. The radiation-pressure force is responsible for the jet acceleration. Right panel: The radial profiles of the gravitational 
force (black), the centrifugal force (magenta), the gas-pressure force (green), the radiation-pressure force (red), and the Lorentz force (blue), at z = 40 
r B . The solid lines and the dashed lines indicate the inward force and the outward force. We understand that the Lorentz force is responsible for the 
collimation of the jet. In both panels each quantity is time-averaged over 1 s, and the shadowed areas indicate the jet region in which the outflow velocity 
exceeds the escape velocity. 
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of the Lorentz force in each region, we display in the right wall 
surface of Figure 1 the color contours of the magnetic pressure 
from the azimuthal component of the magnetic field, B^/8ir, 
overlaid with the Lorentz force vectors (arrows), J x B/pc, of 
the RMHD flow on the meridional plane. 

As we already mentioned, the radial, inward component of 
the Lorentz force, i.e., pinch effect of the magnetic field, dom- 
inates over other forces, including the radiation-pressure force, 
in the jet region. This seems, at first glance, odd, since the 
magnetic energy density is much less than the radiation en- 
ergy density; say, by a factor of ~ 20 at (r, z) = (18 r s , 40 r s ). 
Nevertheless, magnetic force can dominate over the radiation 
force, since the radiation-pressure force is largely attenuated in 
the horizontal direction (in the disk plane) because of highly 
anisotropic distribution of radiation energy. In the high lumi- 
nosity flow with large matter column density, the accumulation 
of material around the black hole is responsible for generating 
rather flat distribution of the radiation energy density in the hor- 
izontal direction. As a result, the horizontal component of the 
radiation-pressure force, which is proportional to the horizon- 
tal gradient of the radiation energy density, is largely reduced. 
By contrast, magnetic energy density shows high concentration 
("magnetic tower") in the cylindrical region surrounding the 
black hole, asserting large inward magnetic force at the inner 
surface of the magnetic tower. The magnitude of the Lorentz 
force can be by ~ 5 times as large as that of the radiation force 
there. 

Note that a magnetic field is by nature intrinsically expan- 
sive. Thus, there must be something to prevent the field from 
expanding sideways, thereby decollimating the whole config- 
uration (Uzdensky & MacFadyen 2006; Spruit 2010). In our 
simulations, it is the geometrically thick accretion flow sup- 
ported by strong radiation pressure force that prevents the filed 
from expansion. We can prove this by looking at Figure 2 
(right) that the central region (the magnetic tower) is confined 
by the combination of the radiation pressure force (the red 
curve) and the gravitational force (the black line). Cautions 
should be taken for the bird's eye view in Figure 1; it only 
traces very high density regions ("core" of the accretion flow) 
and substantial amount of material exists around the magnetic 
tower, as is illustrated in the left wall in Figure 1 . 

We started the simulations from the particular initial con- 
ditions, and we admit that our results may depend on them. 
We need further simulation studies starting from different ini- 
tial magnetic field configurations(Igumenshchev et al. 2003; 
Beckwith et al. 2009). Cases of continuous addition of a mag- 
netic field in a flow should also be necessary to examined 
(McKinney & Narayan 2007). We are planning such simu- 
lations in future work. 

4. Discussion 

Recently, Fender et al. (2004) compiled the X-ray observa- 
tional data of Galactic microquasars in the luminosity-hardness 
ratio diagram and reported a evolutionary paths in the diagram. 
Importantly, they found a common condition for the occurrence 
of relativistic jets in microquasars; that is, relativistic jets only 
appear when the source exhibits a hard-to-soft transition at high 
luminosities. The typical photon luminosities (L p h), at which 



relativistic jets are observed, are close to (or, at least, a few ten 
percents of) Le- Our results are consistent with such observa- 
tional facts. 

The present RMHD simulation, however, does not produce 
highly relativistic jet with the bulk Lorentz factor of a few. 
This is partly because our simulations are restricted to the 
Newtonian dynamics. In future relativistic RMHD simulations 
will able to show really relativistic jets by extracting energy 
from the black hole spin energy. Recently made general rel- 
ativistic MHD simulations can indeed achieve high Lorentz 
factor of jets (Blandford & Znajek 1977; Koide et al. 2002; 
McKinney 2006), though radiation processes are not consid- 
ered in their simulations. In this sense, our simulations are 
more relevant to mildly relativistic, persistent jets from a pecu- 
liar source, SS433. In fact supercritical accretion is indicated 
in this source by number of observations (Margon et al. 1979; 
Kotani et al. 1996; Fabrika 2004). 

To conclude, complex interaction between matter, radiation, 
and magnetic field would be an important factor for the rela- 
tivistic jet formation at high luminosities. Although radiation 
energy dominates, magnetic fields manifest themselves in the 
form of magnetic tower structure, which provides a collimation 
mechanism, which radiation field cannot give. Such situations 
are expected to always occur at above 10% of the Eddington 
luminosity, above which the radiation energy density exceeds, 
or is at least comparable to, internal gas energy density, and 
thus exceeds the magnetic energy density. Our global RMHD 
simulations demonstrate the universal role of magnetic field, as 
well as the importance of radiation-hydrodynamics in physics 
of black hole accretion disks and flows. 
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